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The Mars Global Surveyor spacecraft has been located in a near-circular, polar, and 
low-altitude mapping orbit about Mars for six years, since February 1999. The spacecraft 
is tracked routinely by the antennae of the Deep Space Network (DSN), using the X Band 
radio system of the spacecraft. These tracking data have been used for routine spacecraft 
navigation, and for radio science studies, such as the estimation of the static and time- 
varying gravity field of Mars. In this paper we describe the methodology for reduction of 
these data in order to estimate the Mars atmospheric density (normalized to an altitude 
380 km) over half a solar cycle, where we discern the correlation of the density with the 
incident solar flux, and the 27-day solar rotation. The results show that the density at the 
MGS altitude varies from a mean of 0.7 x 10 -17 grams/cm 3 near aphelion to a mean of 3.0 
x 10 -17 grams/cm 3 near perihelion. 
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Subscripts 

i macromodel plate index 

l spherical harmonic degree 

to spherical harmonic order 

Conventions 

Ctm,Sim Normalized spherical harmonic coefficients 
Pi m Normalized associated Legendre function 

fiCim, 6 Sim Variations in the spherical harmonic coefficients, Ci m and Si m 
Ay Velocity impulse 

n Plate surface unit normal vector 

s Unit vector to the radiation source 
v Satellite velocity vector relative to the atmosphere 

Ai Area of macromodel plate i 
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a dr Spacecraft acceleration due to atmospheric drag 

a pp Spacecraft acceleration due to radiation pressure 

c Speed of light 

C r Reflectivity coefficient 

Fd Drag scale factor 

GM Gravitational constant, m 3 /s 2 

&2 Love number of degree two 

L s Areocentric solar longitude 

m sc Spacecraft mass 

r radius 

R e Reference radius for gravity model 
t Time 

to Reference epoch, January 1.0 2000 
U Gravity potential 
AMD Angular momentum desaturation 
AU Astronomical unit 
DSN Deep Space Network 

ECMWF European Centre for Medium-Range Weather Forecasts 
F10.7 Solar flux at wavelength of 10.7 cm. 

HGA High Gain Antenna 

LMST Local Mean Solar Time 

MGS Mars Global Surveyor 

MJD Modified Julian Date 

MOLA Mars Orbiter Laser Altimeter 

sfu Solar Flux Unit, 10~ 22 W/(m 2 Hz) 

USO Ultra Stable Oscillator 

X MGS spacecraft body-fixed axis 

Y MGS spacecraft body-fixed axis 

Z MGS spacecraft body-fixed axis 

Symbols 

P Specular reflectivity 

6 Diffuse reflectivity 

rii Angle between plate normal, hi, and the satellite velocity vector, v 
A longitude 

ijJi,lo 2 Annual and semiannual Mars orbit periods 
$ Radiation flux at the satellite 
4> latitude 

p Atmospheric density 

cr Standard deviation 

9i Angle between plate normal, hi, and radiation source direction, s 

I. Introduction 

The Mars Global Surveyor spacecraft, the successor to the Mars Observer spacecraft, arrived at Mars in 
September 1997. It heralded the beginning of a new age of Mars exploration, following a hiatus of nearly 16 
years, since the end of the Viking Mission, and nine years since the Soviet Phobos II mission to Mars. Mars 
Global Surveyor had two attributes of prime importance for Radio Science: (1) the spacecraft carried an X 
Band communication system, compared to the S Band tracking system of the Viking Orbiters and Mariner 
9; and (2) the spacecraft was located in a near-circular, and polar orbit at a mean altitude of 400 km. The 
higher frequency of the X Band system (near 8 Ghz) permitted the routine collection of tracking data with 
a noise of 0.10 mm/s (or better), compared to the noisier S Band systems (near 2 Ghz) available in the late 
1970’s. 1 In addition, the near-circular orbit permitted a near-uniform and routine synoptic mapping of the 
planet that was not possible from the Vikings and Mariner 9 since those spacecraft were located in highly 
elliptical near-12-hr and near-24-hr orbits. 

The highly precise data from MGS, and Mars Odyssey, have permitted the development of refined models 
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of the Mars gravity field. 2-5 Solutions have also been developed for the time variations in the gravity field 
of Mars 5-7 and to refine the ephemeris of Mars, 5 ' 8 and to estimate the Mars &2 Love number. 5 ’ 7 ’ 9 

In the reduction of tracking data, it is necessary to model all the forces acting on the orbiting spacecraft. 
Atmospheric drag is modelled using the modified Stewart atmosphere model 10 in the case of GMM2B and 
MGM1041c gravity solutions, 2,4 and a version of the MARSGRAM model 11 in the case of the JPL solutions. 5 
The solved-for drag scale coefficients describe corrections to the underlying atmospheric density model used 
in the precision orbit determination. In this paper we review the derivation of the Mars atmospheric density 
from analysis of MGS data for six years (3 Mars years) since the entry of the spacecraft into the low altitude 
mapping orbit in February 1999. 


II. Data and Methodology 


II. A. Data Description 
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The primary data for Radio Science are the track- 
ing data of Mars Global Surveyor obtained by the 
Deep Space Network (DSN). The observables in- 
clude Doppler and Range data from the three Deep 
Space Complexes at Canberra (Australia), Madrid 
(Spain) and Goldstone (California). 1 The Doppler 
data include both 1-way, 2- way, and 3- way Doppler. 

The 2-way Doppler data involve the same DSN 
station for transmission and reception. The 3- 
way Doppler data involve different DSN stations 
for transmission and reception. The 1-way Doppler 
data originate on the MGS Ultra Stable Oscillator 
(USO), whose primary purpose is to support the 
MGS radio occultation experiment. 1 No 1-way data 
for orbit determination purposes were collected prior 

to March 2000, however afterwards, these data constitute on average 45 percent of the data in a given data 
processing arc. We processed the tracking data in arcs (segments) with an average length of 4.65 days. The 
arclengths were meticulously chosen based on the distribution of the 2-way tracking data, with every arc 
always starting and stopping with a pass of 2-way Doppler data. In addition, a further constraint was to 
select start and stop times so as to avoid AMD events close to the arc endpoints which would be poorly 
determined or aliased with the initial state vector of the arc. 

We included data from the Gravity Calibration Or- 
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Figure 1. Hours of daily Doppler tracking for Mars Global 
Surveyor. 


bit (February 1999), the period of fixed high gain 
antenna mapping (March 1999), and regular map- 
ping mission (after April 1999) through July 7, 2005. 
A total of 401 data arcs were selected for this study. 
We avoided data in the vicinity of solar conjunc- 
tion (July 2000, August 2002, and September 2004) 
since the high noise due to the solar plasma caused 
the quality of the orbit solutions to deteriorate. We 
also did not use tracking data taken while the space- 
craft was in contingency mode or safe mode. The 
quantity of DSN data ranged from near continuous 
tracking during some mission phases (except for oc- 
cultations by Mars as seen from Earth) to as little as 
5 hrs of contact time per day (See Fig. 1). The orbit 
geometry experienced by MGS is depicted in Fig. 2. 
The orientation of the MGS orbit Earth beta angle 
varied from edge-on (orbit plane parallel to the line 
of sight, and Earth beta angle of 0°) to nearly face-on (Earth beta angle of 90° at different phases of the 
mission. The orbit orientation determines the sensitivity of the Doppler data and the recovery of geophysical 
and orbit parameters. For example, the Mars gravitational constant, GM, is best determined during an 



Figure 2. Earth beta angle and solar beta angle for Mars 
Global Surveyor from 1999 to 2005. An Earth beta angle of 
0° corresponds to edge-on orbit geometry, whereas an Earth 
beta angle of 90° corresponds to face-on orbit geometry. 
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edge-on orbit orientation. 2 MGS is located in a sun-synchronous orbit at approximately 2 PM. However, as 
illustrated in Fig. 2, the actual LMST varies by ± 40-50 minutes every Mars year due to the eccentricity of 
the Mars orbit. 

We processed the MGS tracking data using the NASA GSFC GEODYN orbit and geodetic parameter 
estimation program, which is a batch least-squares orbit determination filter that can process both terrestrial 
and interplanetary tracking data. 12 

II. B. Force Modelling 

The reduction of the tracking data requires that we model all the forces that act on the spacecraft orbit. 
We model the Mars gravity field to 90x90 in spherical harmonics using MGM1041c, a gravity model based 
on MGS data through May 2002. 4 * We model the Mars gravity field with normalized coefficients (Cj m , Si m ) 
using the equation 13 


r GM GM v — / Re\ - . . . r=. . .. = . . ..i , . 

U = —y + —y ~ j Plm ( Sln 4) [ c im cos(mA) + Sim sin(mA)] (1) 

where GM is the universal constant of gravitation times the mass of Mars, l is the degree, m is the order, 
Pim are the fully normalized associated Legendre polynomials, R e is the reference radius of Mars (3397 km 
for MGM1041c), (j> is the latitude, and A is the longitude. By definition, the degree one terms are zero, since 
we choose the origin of the coordinate system to be at the center of mass of the planet. 

Other forces modelled include the third body perturbations due to the Sun, planets, and satellites of 
Mars, Phobos and Deimos, and the Mars solid tide using an a priori Love number of 0.10. 

MGS uses momentum wheels to maintain its nominal attitude. Since thruster firings from momentum 
wheel desaturations can perturb the orbit, we adjust accelerations radial, along-track, and cross-track to 
the orbit at the times of the AMD maneuvers. Early in the mapping mission, these maneuvers occurred as 
frequently as 3 to 4 times per day. However, after August 16, 2001, the spacecraft was no longer maintained 
in a strictly nadir orientation (+Z axis with the instrument deck facing nadir). Rather, the Z-axis was pitched 
16° from nadir in the orbit plane. This ’Relay- 16’ configuration had the effect of reducing the frequency of 
the AMD desaturations to typically one per day. 14 

The solar radiation pressure, planetary radiation pressure (due to the planet’s albedo and thermal emis- 
sion), and the atmospheric drag are modelled using a macromodel. 15 The spacecraft is modelled as a series 
of elemental flat plates, and the contributions of each plate are vectorially summed to obtain the total 
acceleration. For the solar or planetary radiation pressure, the total acceleration is given by 


OflP 


y Aj COS 9j 

c r ■G-' m sc 
1=1 



+ 2 Pi cos 9i 


hi + (1 - Pi)s 


(2) 


where $ is the radiation flux at the satellite, c is the speed of light, A, is the area of the 'tth surface plate, n* 
is the surface normal unit vector, s is the unit vector along the direction from the satellite to the radiation 
source (either the Sun or a surface element of the planet), 0* is the angle between the plate normal, hi, and 
the radiation source direction, s, Pi is the specular reflectivity (percent of the total incident radiation), d, is 
the diffuse reflectivity (percent of incident radiation), C r is the reflectivity coefficient (which is estimated), 
and Lambert’s cosine law of diffuse reflection is assumed. 16 

The drag acceleration is modelled as, 


1 

a DR = ~- 


x pv\v\ 


Fd_ 

m sc 


Imax 

Y Ai COS T]i 
i=l 


( 3 ) 


where p is the atmosphere density, v is the satellite velocity vector, |u| is the satellite velocity magnitude, Fd 
is the drag scale coefficient (which is estimated), m sc is the spacecraft mass, Ai is the area of the ith plate, 
and r)i is the angle between the surface normal unit vector, hi, and the velocity vector, v. The spacecraft 
mass is provided by the spacecraft team at Lockheed Martin, and varies from 735.9 kg in February 1999 to 
697.1 kg in June 2005 due to fuel consumption from orbit trim maneuvers and AMD events. 

We derive the a priori atmospheric densities from the modified Stewart atmosphere model. 10 In the 
original Stewart model, the F10.7 solar flux was modelled as an analytic function. Instead, we use the 


4 of 13 


American Institute of Aeronautics and Astronautics Paper 2006-6395 



actual F10.7 flux smoothed over three solar rotations (81 days) to drive the a priori atmospheric density 
model. The planetary radiation pressure uses spherical harmonic models defined over different areocentric 
solar longitudes, L s , to account for the seasonal variations. 17 We use a 12-plate macromodel for MGS: six 
for the spacecraft bus, four for the Y+ and Y- solar arrays (front and back), and two for the high gain 
antenna (front and back). We obtained a priori values from Lockheed Martin, and tuned these using 18 
months of MGS data. 18 In this paper, we describe the results of a new tuning for some of the macromodel 
parameters. The use of a macromodel requires detailed knowledge of the spacecraft attitude and how the 
articulating appendages (the Y+ and Y- solar arrays and the high gain antenna) are oriented with respect 
to the spacecraft and in inertial space. We use spacecraft attitude data in the form of quaternions supplied 
to the MGS project, and we fill in the quaternions during telemetry data gaps. 
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Figure 3. Nonconservative accelerations on MGS, 1999 to 2005. 

In Fig. 3 we depict the nongravitational accelerations calculated from daily averages for MGS from 1999 
through 2005. The mean accelerations are 8.4 x 10 -9 m/s 2 for the planetary radiation pressure, 7.1 x 10~ 8 
m/s 2 for the solar radiation pressure, and 3.8 x 10~ 9 m/s 2 for atmospheric drag. Due to seasonal variations 
and the change in the solar activity over the solar cycle, the atmospheric drag acceleration on MGS varies 
between 6.0 x IIP 10 m/s 2 and 1.0 x 10” 8 m/s 2 , or a factor of 3.3. The radiation pressure perturbations vary 
by a factor of 1.8 between Mars perihelion and aphelion. The amplitude of the accelerations might imply 
that the radiation pressure perturbation will swamp the perturbation due to atmospheric drag. This is not 
the case, as the atmospheric drag acts primarily along-track. In contrast, the planetary radiation pressure 
acts primarily in the radial direction, and for the MGS orbital geometry (sun-synchronous at approximately 
2 PM local mean solar time), the solar radiation pressure acts primarily cross-track and radial to the orbit. 
The mean accelerations for two orbital arcs, one near perihelion, and one near aphelion are summarized in 
Table 1. The ability to estimate a drag scale factor, and hence the density variations will depend on the 
length of period chosen for the estimate of mean density, and the quality of the orbit solutions that are 
obtained. 
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Table 1. Mean Radiation Pressure and Drag Accelerations on MGS for Test Arcs 


Dates 

Perihelion, MJD 51517-51522, 

Aphelion, MJD 52551-52556, 


Along-track 

Cross-track 

Radial 

Along-track 

Cross-track 

Radial 

Solar radiation pressure* 

0.8 

-29.2 

-21.8 

0.8 

-21.8 

-15.8 

Planetary radiation pressure* 

-0.4 

1.8 

10.5 

0.1 

0.9 

7.0 

Atmospheric drag* 

-6.9 

0.1 

0.01 

-2.3 

-0.003 

-0.0007 


* Accelerations are in units of nanometers/s 2 . 


II. C. Measurement Modelling 

The measurement modelling used the DE410 plan- 
etary ephemeris, which incorporated MGS and 
Odyssey data to improve the Mars ephemeris. 8 The 
latest set of DSN station coordinates were used, 21 
where we note that different coordinates must be 
applied to DSN station 65 after March 2005, when 
the station was relocated. The measurement model 
applies geometric Earth tide and ocean loading cor- 
rections. We use in situ weather data collected at 
the three DSN complexes to compute Earth tropo- 
sphere corrections according to a modified Hopfield 
model. 22 Occasionally, these weather data must be 
filtered to replace spurious information. For exam- 
ple, we show in Fig. 4 that the 1999 pressure mea- 
surements from Tidbinbilla were replaced with pres- 
sure data from the ECMWF weather model. 

The relativity modelling includes the relativistic 
delay for both the Sun, Jupiter and Saturn. The 
two-way relativistic delay on Earth-Mars radio sig- 
nals ranges from 4-37 km for the Sun, 0.9 to 2.7 m 
for Jupiter, and 0.1 to 0.5 m for Saturn. 17 In addition, the relativity modelling includes transformations 
between coordinate and atomic time. 19 

We use the Mars IAU2000 reference system. 20 A new Mars reference system has been proposed that uses 
a different sequence of rotation angles and whose parameters are derived from Viking Lander, Pathfinder 
and other data. 5 Work is underway to implement this model in our GEODYN orbit determination program, 
however it was not available to support this present analysis of MGS data. 

II. D. Initial Orbit Determination Results 

The adjusted parameters in a data arc included the spacecraft state, accelerations to model the AMD 
events, a solar radiation reflectivity coefficient ( C r ), daily drag scale coefficients ( Fd ) in most arcs, pass-by- 
pass frequency biases for the 1-way and 3-way data, and pass-by-pass biases for the range data. The average 
RMS of fit was 0.216 mm/s for the 2-way and 3-way data (computed over 401 arcs), and 0.204 mm/s for 
the 1-way data (computed over 326 arcs). Arcs close to opposition had the lowest RMS of fit, 0.14 to 0.18 
mm/s. Even though we did not process arcs within several weeks of solar conjunction, the increased noise is 
still evident in the data (see Fig. 5). The mean reflectivity coefficient, C r was 1.041 with a a of 0.093, but 
exhibited distinct Mars seasonal variations with a minimum of 0.78 and a maximum of 1.26. From Table 1, 
we see that near aphelion, the mean solar radiation pressure acceleration per arc can amount to 30 percent 
of the mean drag acceleration. Thus, if the C r variations are caused by orbit mismodelling, then the orbit 
errors could alias into estimates for the drag scale factors, and hence the estimates of atmospheric density. In 
order to mitigate this potential source of error, we used the MGS data to estimate macromodel reflectivity 
parameters, the Mars &2 Love number, as well as annual and semiannual variations in the leading zonal 
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Figure 4. Actual and ECMWF weather model pressure data 
for the Canberra DSN complex (Tidbinbilla) in 1999. The 
ECMWF data were used to replace spurious in situ mea- 
surements. 
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Figure 5. Doppler RMS of fit for MGS arcs from February 1999 through June 2005. 


harmonics, C 20 and C 30 . We then recomputed the MGS data arcs with these updated parameters. These 
parameter estimations are described in the next section. 

II.E. Estimation of Geophysical and Macromodel Parameters 

A major attribute of the Mars atmosphere is the annual condensation and sublimation of a sizable fraction 
of the atmospheric mass. This process produces variations in the Mars gravity field which can be sensed by 
orbiting spacecraft such as MGS and Mars Odyssey. 5-7 The most significant of these gravity field variations 
are in the zonal harmonics. We choose to model these zonal variations according to the equations 

8 C 2 q = A\ coswi(f — to) + B 1 sinwi(t — to) + A 2 cosw 2 (t — to) + B 2 sinw 2 (t — to) (4) 

and 

<5C3o = D\ coswi(t — to) + Ei sinwi(t — to) + D 2 cos w 2 (t — to) + E 2 sinw 2 (t — to) (5) 

where 8 C 2 q and 8 C 30 are the variations in the normalized zonal harmonics, uj\ and u> 2 are the annual and 
semiannual Mars periods, A\ , B\, A 2 and B 2 are the cosine and sine coefficient components of the annual 
and semiannual variation for SC 2 o- Similarly, D\, E\, D 2 , E 2 are the cosine and sine coefficient components 
of the annual and semiannual variations for 8 C 30 . The reference epoch to is January 1.0, 2000, and the 
argument t — to is expressed in days since the reference epoch. We express the variations with respect to the 
mean values for C 2 q and C 30 from the MGM1041c gravity model. 4 

We estimated numerous combinations of macromodel and geophysical parameters. For the macromodel, 
we only chose those parameters to which the MGS data seemed most sensitive, and we rejected those solutions 
where the macromodel reflectivity estimates were unrealistic. In addition to the periodic zonal terms, 8 C 2 o 
and SC 30 , the Mars k 2 Love number, we also adjusted the GM’s of Mars and Phobos, although the GM 
of Phobos changed negligibly from the a priori value of 720,000 m 3 /s 2 . The MGS data did not have any 
sensitivity to the Deimos GM. We tested solutions for the annual and semiannual variations in 5C 2 1 and 
5 S 2 1 , but found those to be at the edge of significance. The final geophysical parameter estimates are listed 
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in Table 2. In Table 2, we also list the formal solution standard deviation scaled by a factor of ten, and the 
percent variation in each parameter estimate from the different trial solutions. The MGS data determine 
most robustly the SC 30 variations followed by the SC 20 semiannual variation. As long as we use the entire 
set of data (3 Mars years), and we vary only the sets of parameters estimated, with the exception of the < 5(720 
annual terms, the geophysical parameter estimates appear quite stable. If the solutions are attempted using 
two 3-year (approximately 1.4 Mars years) subsets (April 1999 to June 2002; July 2002 to June 2005), the 
recovered values show little consistency with the solutions using the full six years of data. The full six year 
data set is necessary to resolve the estimated parameters, since they all (some by definition) have annual 
or semiannual signals, especially in a sun-synchronous orbit. We note that our recovery of Mars k 2 is quite 
similar to that of Ref. 5. The extremely poor determination of the 5C 20 annual variation is likely due to the 
deleterious effect of the AMD’s which contaminate the parameter recovery. 


Table 2. Estimated Geophysical Parameters from MGS Data, 1999-2005 


Parameter 

A priori value 

Estimate 

Error* 

Percent variation^ 

Ai,< 5 C 2 o annual cosine 

0 

6.218 x 10- 11 

1.3 x 10 - 8 9 

224 

B \ , < 5(720 annual sine 

0 

-4.317 x 10 - 10 

1.5 x 10 " 9 

23 

A 2 ,SC 20 semiannual cosine 

0 

1.531 x 10 ~ 9 

1.4 x 10 " 9 

9.6 

B 2 , SC 2 o semiannual sine 

0 

1.105 x lO " 9 

1.4 x 10 " 9 

14.1 

D \ , JC 30 annual cosine 

0 

2.233 x 10 -9 

8.5 x lO -11 

0.5 

Ei , <56 30 annual sine 

0 

1.946 x 10 “ 9 

7.1 x 10“ u 

0.9 

D 2 , <5630 semiannual cosine 

0 

4.088 x 10 - 10 

7.5 x 10 " 11 

1.2 

E 2 ,SC 3 o semiannual sine 

0 

3.111 x lO " 10 

7.7 x lO ” 11 

5.5 

k 2 

0.10 

0.176 

0.041 

3.9 

Mars GM (m 3 /s 2 ) 

4.282837024 x 10 13 

4.282836798 x 10 13 

1.37 x 10 6 

1.4 x 10-® 


* We report the formal errors of the least squares solution scaled by a factor of ten. 

i The percent variation is derived from examination of the variation in the parameter estimates from multiple solutions. 


Table 3. MGS Macromodel Adjustments 


Plate 

Area 

(m 2 ) 

A priori* 
Reflectivity 

(P, 5) 

Adjusted 

Reflectivity 

03 . «)* 

TX Bus 

3.30 

0.13, 0.52 

5 = 0.3165 

-X Bus 

3.30 

0.13, 0.52 

6 = 0.4008 

+Y Bus 

3.56 

0.13, 0.52 


-Y Bus 

3.56 

0.13, 0.52 


TZ Bus 

2.31 

0.13, 0.52 


-Z Bus 

2.31 

0.13, 0.52 


Solar arrays (front)'*’ 

15.70 

0.0494, 0.1976 

P = 0.0653 

Solar arrays (back)"* 

15.70 

0.0790, 0.2822 

5 = 0.4494 

HGA (front) 

1.94 

0.10, 0.40 


HGA (back) 

1.94 

0.13, 0.52 



* A priori macromodel values from Ref. 18. 

t The Y+ and Y- solar arrays are treated as a unit in the macromodel 
adjustment. 

t Only the adjusted reflectivity values are listed. 


The macromodel reflectivity parameter estimates are given in Table 3. In the MGS spacecraft coordinate 
system during normal mapping operations (before the start of ’Relay-16’ mode), +Z refers to the instrument 
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deck portion of the spacecraft bus and faces nadir, +Y and -Y face the cross-track directions, +X and -X face 
in the velocity and anti-velocity directions. The solar arrays extend along the +Y and -Y axes and normally 
track the Sun. The high gain antenna is located on a boom extending in the -Z direction, and tracks the 
Earth using a dual-axis gimbal. Near the subsolar point, the front sides of the solar arrays face the Sun, and 
the back sides of the solar array face the planet when the planetary radiation pressure perturbation is at 
its peak. The solar array reflectivity values are composites that represent the contribution of the inner and 
outer solar arrays, the drag flaps, and the yoke. 23 The Y+ and Y- solar array parameters are tied together 
since we cannot obtain separate estimates for each panel. Prom previous experience, we know that the most 
important parameter to adjust in order to reduce macromodel error is the specular reflectivity of the solar 
array. 24,25 The diffuse reflectivity of the back side of the solar arrays was adjusted to accommodate errors 
in the planetary radiation pressure model. An analysis of the information matrix prior to inversion also 
indicated that the +X and -X panels had heightened sensitivity, so the diffuse reflectivities for those plates 
were also adjusted. 


II. F. Results of Updated Orbit Computations 



In Table 4, we give the results of two additional 
series of orbit computations for the MGS data arcs: 

(1) applying the new time- variable gravity estimates 
and the new value of the k 2 Love number, and (2) 
applying both the new geophysical parameter esti- 
mates and the new macromodel reflectivities. The 
application of the SC 2 o and SC 30 variations and the 
new k 2 improve the average RMS of fit by about 
10 percent or 0.02 mm/s. The subsequent applica- 
tion of the macromodel has minimal impact on the 
average level of fit to the data. However, the most 
notable result is the reduction in the scatter of the 
reflectivity coefficient, C r , as measured by the stan- 
dard deviation computed over all 401 data arcs. The 
standard deviation, a is reduced from 0.093 in the 
original processing, to 0.065 with the application of 
the new geophysical parameters, and to 0.050 with 
the application of both the geophysical and macro- 
model parameters. We show the reflectivity coef- 
ficient in Fig. 6 for the original and the new time 
series. The strong annual excursions in the reflec- 
tivity coefficient, particularly near MJD 52700, have been reduced. The average RMS difference in the orbits 
between the original and updated processing is 0.17 m in the radial direction, 4.17 m cross-track, and 2.63 
m along-track. Sometimes the orbit differences between the new and updated processing can be large. For 
example, the cross-track orbit differences are 20-50 m whenever the Earth beta angle for the MGS orbit 
is less than 20°. This is understable as during an edge-on orbit orientation (Earth beta angle of 0°), the 
Doppler data have no cross-track sensitivity, and changes in the force model cross-track to the orbit will 
have maximum effect. 
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Figure 6 . Radiation pressure reflectivity coefficient, C r » for 
the original and the updated time series of MGS data arcs 
from 1999 to 2005. The updated time series includes the an- 
nual and semiannual variations in the zonal harmonics, 6 C 20 
and SC 301 the new k% estimate, and the updated macromodel 
reflectivity parameters. 


III. Atmospheric Density Results 

We first review the drag scale factors obtained by GEODYN using the final updated orbit time series 
(TVG + k 2 + macromodel in Table 4). The daily drag scale factors were filtered to remove some with 
negative values, those whose relative errors were greater than 50 percent, or whose standard deviations were 
excessively large. We retained 1760 drag scale factors out of 1794, and show the relative error in percent in 
Fig. 7. The relative error is derived from the formal error of the drag scale factors in each orbit determination 
solution, and reflects the a priori data weighting (discussed earlier) as well as the distribution and quality of 
the tracking data in each data arc. It is evident that the daily drag scale factor estimates are much noisier 
early in the mission, when the spacecraft experienced 3-4 AMD maneuvers per day. Through August 16, 
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Table 4. RMS of Fit Summary for MGS Arcs from 1999 to 2005 

Orbit Series RMS, 2-way & RMS, 1-way Reflectivity 

3-way Doppler Doppler Coefficient 

(mm/s)* (mm/s)* C r 



Mean 

Median 

Mean 

Median 

Mean 

Median 

Std. Dev.* 

Original 

0.214 

0.198 

0.203 

0.186 

1.042 

1.037 

0.093 

+TVG + k 2 * 

0.195 

0.182 

0.180 

0.165 

1.048 

1.040 

0.065 

+TVG + k 2 + new macromodel* 

0.195 

0.182 

0.181 

0.164 

1.059 

1.054 

0.050 


‘The statistics are computed over 401 arcs for the 2-way and 3-way Doppler, and the reflectivity coefficient, C r . The 
statistics for the 1-way Doppler are computed over 326 arcs, 
t Standard deviation. 

tThe second time series, TVG + k 2 , applies the SC 20 and 5 C 30 variations and the new k 2 - The third time series, 
TVG+fc 2 +macromodel, applies both the new geophysical parameters and the new macromodel reflectivity values. 


2001 or the period of standard mapping, the mean relative error is 6.9%, the median is 5.4%, and the RMS 
is 9.2%. In contrast, during ’Relay-16’ mode, the mean relative error is 2.4%, the median is 1.7% and the 
RMS is 3.6%. These relative errors are important as they will map directly into the derived densities. 
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Modified Julian Date 

Figure 7. Percent relative error in the drag scale coefficients, Fo, estimated by GEODYN for MGS from 1999 to 2005. 

The recovered densities are shown in Fig. 8 from 1999 through 2005. We use the 10.7 cm radio flux as 
a proxy for the EUV variability that is known to drive the behaviour the planetary thermospheres. This 
flux, which is measured at the Earth, is corrected by scaling the daily F10.7 defined at 1 AU to the actual 
Mars distance from the Sun, and the variations in the Sun-Earth-Mars angle as the planets move about the 
Sun. 26 The period from 1999 to 2005 includes the last maximum of the solar cycle 23, and the approach to 
the current solar minimum. The combination of variations in solar flux, and the distance of Mars from the 
Sun cause the density at the MGS altitude to vary by a factor of five, from 1 to 5 x 10~ 17 grams/cm 3 . The 
mean density near aphelion (L s = 75°) is about 0.7 x 10~ 17 grams/cm 3 whereas near perihelion (L s = 255°) 
the mean density is 3.0 x 10 -17 grams/cm 3 (see Fig. 9). In addition, regional and global dust storms which 
generally occur within a few months of perihelion, also contribute to the density variability depicted in Fig. 8 
and Fig. 9. 

The density data show evidence of the Mars atmosphere responding to the changing solar flux, induced 
by the near 27-day solar rotation. In Fig. 10, we illustrate the behaviour of the Mars thermosphere over eight 
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month time span. The recovered atmospheric density residuals clearly contain the same 27-day periodicity 
that is evident in the corrected FI 0.7 data. 
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Figure 8. Normalized density at 380 km for MGS and corrected daily solar flux (F10.7). A solar flux unit (sfu) has 
units of 10 ~ 22 W/(m 2 Hz). 
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Figure 9. MGS normalized density at 380 km versus areocentric solar longitude, L s • 


IV. Conclusion 

We have computed the mean density at the MGS altitude over 3 Mars years, from 1999-2005. The 
variations in density are controlled by variations in solar flux, and modulated by both the near 27-day solar 
rotation, and change in solar flux over the solar cycle. Mean densities at the MGS altitude can vary by a 
factor of five between perihelion and aphelion. 

The recovery of atmospheric density requires precise modelling of all forces acting on the MGS spacecraft. 
The macromodel that is used to model the variable cross-sectional area and provide for more sophisticated 
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MGS normalized density residuals (380 km) vs corrected daily FI 0.7 residuals 



modified Julian date 1 950 

Figure 10. MGS normalized density residuals (380 km) vs. corrected daily F10.7 residuals for May 23, 2003 to March 
18, 2004. 


drag and radiation pressure modelling, is a major improvement over using a simplistic model such as a 
cannonball. However, the macromodel does not allow for self-shadowing of spacecraft components, nor does 
it account for spacecraft thermal emission or the radiation interaction between different spacecraft surfaces. 
The a priori modelling of the zonal SC20 and SC30 variations appears to be essential in order to achieve the 
best possible orbits, yet in our analysis we did not obtain a reliable solution for the SC20 annual variations. 
If a more parsimonious strategy were followed to model the AMD’s, such as using information supplied 
by the spacecraft team, or limiting the estimation of accelerations or Ay’s to the axes along which the 
AMD occurred, better solutions for the SC 20 variation might be obtained. The addition of MOLA altimeter 
crossovers as an orbit determination data type could improve the orbit determination solutions for MGS 
through June 2001 and consequently would provide more robust estimates of the MGS density early in the 
mapping mission. 

A logical extension of this analyis being pursued by the authors is to incorporate these MGS density 
data, along with improved solar activity and seasonal variations, and dust storm effects, into an empirical 
density and temperature model for the Mars thermosphere, such as DTM-Mars. 18 
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